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ABSTRACT 


Atmospheric  dispersion  is  shown  to  have  a  significant  effect  on  wideband  coherent  laser  radars 
operating  on  a  transition  of  ordinary  carbon  dioxide  (  CO2 ).  Calculations  are  performed  for  a 
hypothetical  ground-based  laser  radar  that  is  observing  targets  in  low  Earth  orbit.  Linear  frequency 
modulated  (  LFM  )  waveforms  of  bandwidths  200  MHz,  500  MHz,  1  GHz,  and  2  GHz  are  used. 
The  use  of  a  laser  radar  operating  with  a  different  carbon  isotope  ( ^^C^^02  )  is  suggested  to  avoid 
the  dispersion  problem,  and  also  to  reduce  the  atmospheric  absorption. 
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1.  INTRODUCTION 


Carbon  dioxide  ( CO2 )  lasers  can  be  used  in  coherent  optical  radars.  Some  COj  radars  will 
operate  through  the  Earth’s  atmosphere. 

Atmospheric  CO2  and  other  gases  will  absorb  the  laser  radiation.  Because  this  is  resonant 
absorption,  dispersion  will  also  occur.  The  combined  effects  of  absorption  and  dispersion  can 
significantly  degrade  the  resolution  of  the  radar.  These  effects  are  calculated  in  this  paper.  The 
particular  example  chosen  for  calculation  is  for  a  radar  that  is  used  to  determine  the  precise  range  to 
a  satellite  in  Earth  orbit.  This  example  was  chosen  to  show  the  effects  of  CO2  absorption  both  at 
low  and  high  altitudes  in  the  atmosphere. 


2.  METHOD  OF  CALCULATION 


The  following  discussion  describes  the  method  used  to  estimate  the  atmospheric  absorption 
and  dispersion  effects.  Some  of  the  approximations  employed  are  discussed. 

The  Air  Force  Geophysics  Laboratory’s  code  for  computing  atmospheric  spectral 
transmittance  with  high  resolution  ( FASCOD2^ ),  with  its  extensive  compilation  of  spectral  data,  is 
a  convenient  tool  for  calculating  atmospheric  absorption.  In  principle,  the  code  could  be  modified 
to  integrate  optical  path  length  as  well  as  absorption  through  the  atmosphere;  an  approach  similar  to 
this  was  employed  by  Halliday^  to  estimate  the  dispersion  effects. 

Alternately,  the  Kramers-Kronig  relations^  can  be  used  to  estimate  the  dispersion  from  a 
calculation  of  only  the  absorption.  This  method  was  selected  because  of  its  ease  of  use.  The  basic 
form  of  the  relation  used  is: 


(1) 


where  x'  is  the  real  part  of  the  medium’s  complex  susceptibility  ( X  =  X'  +  j  X " )  a  function  of 


angular  frequency  co,  and  x"  is  the  imaginary  part.  (  P.  V. )  indicates  Cauchy  principal  value.  If 
Rg  is  the  complex  refractive  index  of  the  medium  ( =  n  +  j  k ), 


(2) 


1 


but  for  the  atmosphere 


(3) 


I  - 1  I  «  1 


so  that 

n^-1  =  271 X 


(4) 


The  relation  between  the  absorption  coefficient  (  a  )  and  the  imaginary  part  (  k  )  of  a 
medium’s  complex  index  (  )  is 

a  =  47C  V  k  (5) 


where  v  is  the  frequency  (  wavenumber )  of  the  radiation  ( co  =  27t  v  ).  For  a  large  absorption 
coefficient  of  1  cm‘^  and  for  v  =  1000  cm'^  (wavelength  10  pm),  k  =  8  x  10'^.  A  more  typical 
absorption  coefficient  in  the  atmosphere  would  be  of  order  1  km  ^  or  1  x  10’^  cm  ^  Under  the 
approximation  in  Equation  3,  Equation  1  may  be  rewritten  as 


n ( v)- 1 


=  -P.V.  fdv' 

7t  J 


k(v’) 

v’  -  V 


(6) 


Define  a  function  k  ( v  ) 

k(v)  =  Jds  k(v)  (7) 

P 

( defined  as  a  line  integral  along  path  p,  which  may  be  curved  in  the  atmosphere  due  to  refraction ) 
and  another  as 

Ti  (v)=Jds  [n(v)-1  ]  (8) 

p 

then  q  may  be  expressed  in  terms  of  k  by 

q  (v)=-P.  V.  fdv’-!^^  (9) 

7t  J  Y*  .  V  '  ' 
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The  function  k(  v  )  is  obtained  from  the  FASCOD2  output.  In  practice,  the  range  in  v  must 
be  finite,  but  a  suitably  large  interval  may  be  chosen  to  obtain  a  good  estimate  of  the  effect  of  the 
atmospheric  dispersion  near  a  spectral  feature. 

Fourier  methods  may  conveniently  be  used  to  evaluate  Equation  9.  If  the  Fourier  transform 
of  K  ( V  )  is  defined 

oo 

K(x)  =  Jdve^^*'^K(v)  ("lO) 


and  the  Fourier  transform  of  T|  ( v  )  is  similarly  defined  to  be  H  (  x  ),  then  Equation  9  can  be 
expressed  in  the  transform  domain  ( c.f.  Bracewell  )  as 

H(x)  =-i-  sgn  (x)  K(x)  (11) 

where  the  signum  function  is  defined 

f  -1  ;  x<0 

sgn (x)  =  i  0  :  x  =  o  (12) 

[  1  :  x>0 

The  Fourier  transforms  are  performed  in  practice  using  a  fast  Fourier  transform  (  FFT  ). 
Some  care  is  required  to  minimize  the  effects  of  the  finite  interval  of  the  calculation. 

As  a  test  case,  a  Lorentzian  line  shape  is  used  for  k  and  the  dispersion  curve  (  n-1  )  is 
computed  using  the  above  method  and  compared  to  the  true  values.  The  functions  for  k  and  (  n-1  ) 
are 

k(v)= -  (13) 

[(v-v„)  )1 


and 
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(14) 


[(v-Vq  )^  +  (^)] 


where  Kq  is  the  on-resonance  value  of  k,  y  is  the  half  width  half  strength  frequency  for  the  line, 
and  Vq  is  the  line  center  frequency.  For  Kq  =  1  =  y,  and  Vq  =  0,  Figure  1  shows  the  values  for  k 
and  (  n-1  ).  To  test  the  validity  of  the  approximation  to  the  Kramers  -  Kronig  calculation,  the 
values  from  Equation  13  were  used  in  a  calculation  and  compared  to  the  analytical  result  in 
Equation  14.  Figure  2  shows  the  arithmetic  difference  between  the  function  from  Equation  14  and 
the  approximation.  For  the  Kramers-Kronig  calculation,  the  spectral  width  of  the  calculation  was 
100  full  widths  (  200  y )  over  1024  points,  with  an  additional  1024  points  appended  to  make  the 
number  of  elements  in  the  fast  Fourier  transform  equal  2048  points.  The  error  relative  to  the 
analytical  result  is  neglible  here,  because  phase  errors  of  less  than  2  mrad  peak  -  to  -  peak  cannot 
significantly  increase  the  range  sidelobes  or  broaden  the  main  lobe. 

3.  KRAMERS-KRONIG  CALCULATIONS  FROM  FASCOD2  DATA 

Figure  3  shows  the  optical  depth  of  the  atmospheric  absoiption  [  47c  v  k  {  v  )  ]  for  a  path  to 
space  as  calculated  using  the  FASCOD2  program.  The  path  was  from  sea  level,  for  the 
mid-latitude  summer  model,  at  a  zenith  angle  of  45°.  The  strong  absorption  feature  is  that  of  the 
I  P(  22  )  line  of  atmospheric  CO2  at  942.3833  cm‘^  (  c.f.  Bradley  et  al  ^  ).  The  FASCOD2 
calculation  has  been  compared  with  measured  data  from  Kitt  Peak^  and  found  to  give  a  good 
estimate  of  the  CO2  absorption.  Figure  4  shows  the  results  of  a  similar  calculation  but  with  the 
path  from  a  10  km  altitude  to  space.  The  10  km  altitude  represents  a  hypothetical  case  where  the 
radar  would  operate  from  an  aircraft  to  reduce  the  atmospheric  absorption.  Figure  3  shows 
pressure  broadening  of  the  CO2  line  at  low  altimdes,  plus  some  broadband  absorption  that  is  due  to 
water  vapor. 

The  effects  of  the  finite  spectral  interval  in  the  absorption  are  minimized  by  extending  the 
spectrum  with  dummy  values.  Figure  5  shows  the  spectrum  of  Figure  3,  divided  by  the  frequency 
( v  )  and  padded  with  the  values  from  the  ends  of  the  FASCOD  calculations.  Note  that  the  FFT 
will  treat  the  interval  as  one  cycle  of  a  periodic  spectrum,  so  the  left  side  padding  wraps  over  to  the 
right  of  the  interval.  The  padding  also  extends  the  number  of  points  to  a  power  of  2  for  most 
efficient  use  of  the  FFT  algorithm. 
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Figure  1.  k  and  (n-1)  for  a  Lorentzian  line. 
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Figure  2.  Difference  between  ( n-1 )  for  a  Lorentzian  line  and  an  approximate 
Kramers-Kronig  calculation  from  k. 
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Figure  3.  Optical  depth  of  the  atmosphere  computed  using  FASCOD2  for  a  path 
to  space  from  sea  level  at  a  45°  zenith  angle . 
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Figure  4.  Atmospheric  optical  depth  calculation  from  a  10  km  altitude  to  space. 
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Figure  5.  Optical  depth  from  sea  level,  with  padding . 
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Padding  the  spectral  interval  gives  the  approximate  Kramers-Kronig  result  shown  in 
Figure  6.  The  disturbance  caused  by  the  discontinuity  in  the  padded  spectrum  is  set  to  be  as  far  as 
possible  from  the  center  of  the  spectral  data  to  minimize  the  errors  it  introduces.  Most  of  the  small 
residual  error  will  be  in  the  form  of  a  linear  phase  error,  corresponding  to  a  small  time  shift.  The 
computed  phase  function  is  then  taken  only  over  the  spectral  interval  of  the  FASCOD  data. 
Figure  7  is  the  result  for  the  optical  depth  shown  in  Figure  3. 

Figure  8  shows  the  optical  path  difference  curves  [  27t  v  T|  (  v  )  ]  computed  with  the  optical 
depth  in  Figure  4. 


Given  the  function  k  (  v  )  and  the  dispersion  T|  ( v ),  the  transmittance  transfer  function  is 
defined 


T(v)=e 


j2nv(Ti(v)  +  iK(v)] 


(15) 


Figure  9  shows  the  amplitude  and  phase  of  the  transfer  function  computed  from  the  data  in 
Figures  3  and  7.  The  zero  frequency  point  is  taken  to  be  the  CO2  line  center.  This  is  for  a  single 
pass  through  the  atmosphere.  For  an  object  in  Earth  orbit  ( Figure  10  ),  the  transfer  function  for 
the  return  path  of  a  radar  pulse  will  be  displaced  by  the  Doppler  shift  from  the  target,  which  can  be 
as  large  as  1.5  GHz  for  10  |im  radiation.  Figures  1 1  and  12  show  the  amplitude  and  phase  of  the 
two-way  transfer  functions  ( referred  to  the  center  frequency  in  the  target  return  )  for  a  1.0  GHz 
Doppler  shift  with  the  atmospheres  used  in  Figures  3  and  4.  The  1  GHz  shift  from  the  target  is 
sufficient  so  that  the  transfer  function  will  be  dominated  by  the  upward  path  for  laser  waveform 
bandwidths  much  less  than  2  GHz. 


4 .  LINEAR  FREQUENCY  MODULATION 

The  type  of  waveform  employed  for  the  calculations  uses  linear  frequency  modulation 
( LFM )  to  obtain  range  resolution.  The  waveform  is  known  as  a  chirp  .  In  the  complex  notation, 
a  single  chiip  is  represented  as 

j27l[Vot+l{  — )t^] 


where  Vq  is  the  carrier  frequency  (  *  3  x  10*^  Hz ),  B  is  the  chirp  bandwidth,  and  At  is  the  chirp 
duration.  For  the  simplest  case, 
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Figure  6.  Padded  optical  path  difference  computed  by  the  Kramers-Kronig 
method  for  a  path  from  sea  level . 
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Figure  7.  Optical  path  difference  from  sea  level . 
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Figure  8.  Optical  path  difference  for  the  path  from  10  km. 
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Figure  9.  One-way  atmospheric  transfer  function  from  sea  level . 
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Figure  10.  Illustration  of  the  effects  of  target  motion  for  a  satellite  in  Earth  orbit. 


Figure  11.  Transfer  function  for  a  two-way  calculation  from  sea  level,  for  a 
1  GHz  target  Doppler  shift. 
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Figure  12.  Transfer  function  for  a  two-way  calculation  from  10  km,  for  a  1  GHz 
target  Doppler  shift. 
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(17) 


a(t)  =  rect(  — ), 
At 


where 


f  1  ;  I  X  I  <0.5 
rect(x)  =  •{ 

I  0  :  I  X  I  >  0.5 


(18) 


The  pulse  may  be  compressed  by  matched  filtering  (  taking  its  cross-correlation  with  a 
reference  chirp  that  is  the  complex  conjugate  of  the  chirp  itself ),  The  squared  modulus  of  the 
correlation  is 

c(t).(AL±LLf  sinc^[(At -I  t  I  )(— )tl  (19) 

At  At 


for  (  I  t  I  <  A  t  ;  zero  elsewhere )  where 

sine  ( X ) 


sin  ( 71 X  ) 

71  X 


(20) 


For  large  time  bandwidth  products  (  B  At  »  1  )  note  that  Equation  19  can  be  well 
approximated  for  small  times  (  I  t  I  «  A  t  )  by 

c(t)  =  sinc^(Bt)  (21) 

The  first  time  (  corresponding  to  target  range  )  sidelobe  for  the  sine  -  squared  function  is 
13.5  dB  (  0.045  )  below  the  central  maximum.  Hamming  weighting  (  c.f.  Rihaezek^)  will  be 
used  to  reduce  the  range  sidelobes  ( with  some  loss  of  resolution )  throughout  the  remainder  of  this 
discussion.  The  form  of  the  weighting  is: 

W{v)  =  rect(-g-2)  { 0.54  +  0.46 cos [ - g— ^1)  (22) 


5.  RANGE  SIDELOBES 


The  target  matched  filter  response  may  be  expressed  as 

p(t)  =  I  F.T/'  {  V)  U^-{v)  W(  V)  T(v)  }l^  (23) 

where  (  F.  T.*^  )  indicates  inverse  Fourier  transform,  is  the  Fourier  transform  of  the 

transmitted  chirp,  Uj^  is  the  transform  of  the  reference  chirp  ( *  indicates  complex  conjugate),  W  is 
the  weighting  function,  and  T  is  the  two-way  atmospheric  transfer  function. 

Figures  13-16  show  the  matched  filter  responses  for  the  atmospheric  absorption  of  Figure  3 
( along  a  path  to  space  from  sea  level )  for  chirp  bandwidths  of  200  Mhz,  500  MHz,  1  Ghz,  and 
2  GHz,  and  a  Doppler  shift  from  the  target  of  1  GHz  (  transfer  function  of  Figure  11). 
Figures  17  -  20  are  similar  curves  given  a  radar  at  10  km  altitude  (  atmosphere  of  Figure  4, 
transfer  function  of  Figure  12  ).  In  Figures  13-20  the  curves  are  each  shown  with  a  reference 
Hamming  weighted  impulse  response  for  comparison. 

6 .  USE  OF  ISOTOPIC  COj  LASERS 

To  avoid  both  the  absorption  and  dispersion  associated  with  atmospheric  COj,  a  radar  can  be 
operated  on  a  transition  of  a  different  isotopic  combination.  Ordinary  is  ,  which  is 

represented  in  brief  notation  ( Reference  5  )  as  626.  The  isotope  combination  636(  ^^C'^02 )  is  rare 
in  the  atmosphere.  Figure  21  shows  data  from  Reference  6  over  a  spectral  region  containing  two  of 
the  strongest  lines  of  this  isotope  combination.  The  spectrum  is  at  a  resolution  of  approximately 
5  X  10*^  cm‘^  and  is  through  less  than  2  air  masses.  Note  that  the  absorption  due  to  naturally 
occuring  ^^€^^©2  is  very  small,  and  that  there  are  no  accidental  overlaps  with  absorption  features 
due  to  other  species. 
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Figure  13.  Matched  filter  response  for  a  linear  FM  chirp  pulse  with  a  200  MHz 
bandwidth,  for  the  two-way  transfer  function  from  sea  level. 
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Figure  15.  Matched  filter  response,  1  GHz  bandwidth,  from  sea  level. 
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Figure  16.  Matched  filter  response,  2  GHz  bandwidth,  from  sea  level. 
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Figure  1 7.  Matched  filter  response  for  a  linear  FM  chirp  pulse  with  a  200  MHz 
bandwidth,  from  10  km  ( two-way  ). 
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Figure  18.  Matched  filter  response,  500  MHz  bandwidth,  from  10  km. 
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Figure  19.  Matched  filter  response,  1  GHz  bandwidth,  from  10  km. 
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Figure  21.  Section  of  solar  absorption  spectrum  showing  transmittance  of  ^^€^^©2 
laser  energy. 
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7.  CONCLUSIONS 


We  have  shown  that  the  atmospheric  dispersion  effects  on  high-resolution  CO2  laser  radars 
are  significant  beyond  the  simple  loss  of  power  due  to  the  associated  absorption.  The  principal 
effect  is  to  increase  the  range  sidelobes  of  the  radar  to  levels  that  may  be  unacceptable  for  some 
applications.  Range  resolution  is  also  degraded,  but  this  effect  is  small.  The  calculation  on  a  path 
from  a  10  km  altitude  shows  that  the  dispersion  effects  cannot  be  reduced  to  low  levels  by 
operating  a  CO2  laser  radar  from  an  aircraft,  although  the  atmospheric  transmittance  will  be  notably 
improved. 

Similar  effects  would  apply  to  any  CO2  laser  radar  operating  with  a  resolution  ( bandwidth ) 
comparable  to  any  of  those  in  the  above  calculations.  These  would  include  systems  that  obtain  high 
range  resolution  by  using  short  pulses  (0.5-5  ns  duration )  or  phase  coding  schemes. 

The  strongest  lasing  transitions  of  ordinary  CO2  appear  to  be  poor  choices  for  operation 
through  the  atmosphere.  Other  transitions  that  would  be  weakly  absorbed,  especially  those  of  rarer 
isotope  combinations  ( e.g.,  ^^C^^02,  ^^C**02  ),  will  be  preferred. 
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